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Abstract 

We present new results for the virial coefficients Bk with k < 10 for hard spheres in 
dimensions D = 2, • • • ,8. 
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1 Introduction 

The low density virial expansion of the pressure 

P 



Y,B k p k with B 1 = 1 (1) 



k B T , , 

for the hard sphere gas of particles of diameter a in D dimensions defined by the two body potential 

"w-{r tit: < 2 > 

is one of the oldest systems studied in statistical mechanics. The problem was first studied an- 
alytically by by van der Waals PP, Boltzmann j^j, and van Laar |2j who computed the coeffi- 
cients up through B4. The computation of B4 for D = 2 was first done in 1964 by Rowlinson 0] 
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and Hemmer and very recently these analytic computations for B4 have been extended to 
D = 4, 6, 8, 10, 12 by the present authors |HJ, and by Lyberg [J for D = 5,7, 9, 11. 

All other computations for the hard sphere gas are by means of computer. This work was 
initiated in the 1950s for hard discs by Metropolis et al. ,8, and for hard spheres by Rosenbluth 
and Rosenbluth Subsequently B§ and B7 were computed by Ree and Hoover [TUI II 11 IT2*] during 
the 1960s and B% was computed by Janse van Rensberg J2j m 1993. Computations for D > 3 were 
initiated in 1964 by Ree and Hoover ^1] who computed B4 for D = 4, • • • , 11. The coefficients B$ 
and Bq for D = 4 and 5 were computed by Bishop, Masters, and Clarke ^H] in 1999, and Bishop, 
Masters, and Vlasov ^1 have recently calculated B7 in D = 4, 5 and Bg in D = 4. 

In a series of papers ^| EI DD] we have extended these numerical computations by computing 
virial coefficients up through Bio in -D = 2, 3, • • • ,8. We use the method of Ree-Hoover diagrams 
as evaluated by Monte Carlo integration. The details are given in Our results are given in 
TableHJin the form of Bk/B^ 1 where 

It is well known that for hard spheres in D dimensions that for some sufficiently large k which 
depends on D that some Ree-Hoover diagrams for Bk vanish identically for geometric reasons. We 
give (a lower bound on) the number of non- vanishing Ree-Hoover diagrams in Table |21 



Table 2: Number of Mayer and Ree-Hoover integrals 



Order 





2 


3 


4 


5 


6 


7 


8 


9 


10 


Mayer 


1 


1 


3 


10 


56 


468 


7123 


194066 


9743542 


RH 


1 


1 


2 


5 


23 


171 


2606 


81564 


4980756 


RH/Mayer 


1 


1 


0.667 


0.500 


0.410 


0.365 


0.366 


0.420 


0.511 


RH, D — 1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


RH, D — 2 


1 


1 


2 


4 


15 


73 


^647 


>8417 


>110529 


RH, D — 3 


1 


1 


2 


5 


22 


161 


>2334 


>60902 




RH, D = A 


1 


1 


2 


5 


23 


169 


>2556 


>76318 





2 Behavior of Bk for large k 

It is apparent from Tabled that negative virial coefficients occur. This was first observed for B^ 
in ^2]. We observe that because B4 changes sign between D = 7 and D = 8, Bq changes sign 
between D = 5 and 6, and B$ and Big change sign between D = 4 and D = 5. This suggests that 
for large k the coefficient Bk may become negative for dimensions smaller than 5. In particular if 
for D — 2 or D = 3 there were a value of k such that Bk changed sign then approximate equations 
of state obtained from the first ten virial coefficients would be wholly inadequate to obtain the 
radius of convergence of the virial series. 

The most important property of the virial coefficients Bk is not their actual numerical values 
for k less than some finite number but rather their asymptotic behavior as k — > 00 because it is 
the asymptotic value which determines the radius of convergence. Of course no finite number of 
virial coefficients can give information on the k — > 00 behavior unless there is some a priori reason 
to expect that the values of k are already in the asymptotic k — > 00 regime. 
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D 


B 3 /Bl 


B4/B2 


B5/B2 


B e /B 2 




Bs/B 7 2 


B9/B2 


Bio/ B® 


2 


0.782004- ■ ■ 


0.53223180 ■ ■ ■ 


0.33355604(1)* 


0.1988425(42) 


0.1148728(43) 


0.0649930(34) 


0.0362193(35) 


0.0199537(80) 


3 


0.625 


0.2869495 • ■ ■ 


0.110252(1)* 


0.03888198(91) 


0.01302354(91) 


0.0041832(11) 


0.0013094(13) 


0.0004035(15) 


4 


0.506340 • ■ ■ 


0.15184606 ■ ■ ■ 


0.0357041(17) 


0.0077359(16) 


0.0014303(19) 


0.0002888(18) 


0.0000441(22) 


0.0000113(31) 


5 


0.414063 ■ ■ ■ 


0.0759724807 ■ ■ ■ 


0.0129551(13) 


0.0009815(14) 


0.0004162(19) 


-0.0001120(20) 


0.0000747(26) 


-0.0000492(48) 


6 


0.340941 • ■ ■ 


0.03336314 ■ ■ ■ 


0.0075231(11) 


-0.0017385(13) 


0.0013066(18) 


-0.0008950(30) 


0.0006673(45) 


-0.000525(16) 


7 


0.282227 • ■ ■ 


0.00986494662 ■ ■ • 


0.0070724(10) 


-0.0035121(11) 


0.0025386(16) 


-0.0019937(28) 


0.0016869(41) 


-0.001514(14) 


8 


0.234614 ■ ■ ■ 


-0.00255768 ■ ■ ■ 


0.00743092(93) 


-0.0045164(11) 


0.0034149(15) 


-0.0028624(26) 


0.0025969(38) 


-0.002511(13) 



Table 1: Numerical values of virial coefficients. Values for B7 D > 5, B% D > 4, Bg, and Bio are new, and other values improve on 
published literature results for B 5 and higher except for the results for B 5 for D = 2, 3 which are due to Kratky 20 . 



We see in Table the dramatic effect that the number of non-zero Ree- Hoover integrals in two 
dimensions is far less than that of the number of biconnected graphs with non-zero star content. At 
k = 10 in D = 2 we estimate that only 0.022 of the Ree-Hoover diagrams with non-zero coefficients 
have non-zero integrals. 

The dramatic (at least in D = 2) reduction as k — > oo in the number of non-vanishing Ree- 
Hoover diagrams gives a criteria for the size of k needed for to be in the asymptotic region. 

Criteria 1 

The number of nonzero Ree Hoover diagrams has approached its large k behavior. 

For k = 10 this criteria may only be fulfilled for D = 2 and is surely not fulfilled at all for 
D > 5. 

Our second criteria has been presented in our previous paper |17j 
Criteria 2 

The loose packed diagrams (defined to be those with the number of / bonds near their maximum 
value) numerically dominate Bk as k — > oo. 

The validity of this criteria has been studied in detail in ^7]. Here it was seen that for D = 3 
and k > 12 the criteria is well satisfied and that as D increases the criteria is satisfied for smaller 
values of k. However, for D = 2 the criteria was not satisfied even for k as large as 17. 

We thus conclude that there is no dimension in which both of these criteria are simultaneously 
fulfilled though in D = 3 and D = 4 it is possible that they both could hold for some moderate 
values of k such as 12 — 14. 

3 Ratio Analysis 

Even though we have argued that k = 10 may not be sufficiently large to see the true asymptotic 
behavior of Bk it is still of interest to determine what results are obtained if well known methods 
are used to estimate the radius of convergence from the first ten virial coefficients. 

One such way of estimating the radius of convergence is the analysis of the ratios of coefficients 
[211221 where we plot Bkp cp / Bk-i versus 1/fc (and we have normalized the virial coefficients to 
the density p cp of the closest packed lattice). The ratio extrapolated to 1/fc — > will give the 
radius of convergence of the series pa which may also expressed in terms of the packing fraction 
T) = B2p/1 D ~ 1 ■ If the slope of the interpolated points approaches zero for large k then the leading 
singularity is a pole on the positive real axis, if the slope is non-zero then the divergence is algebraic. 

The plot of BkPcp/ 'Bk-i versus ^ for D = 2 is given in Fig.^ Here we observe smoothly falling 
ratios which extrapolate to a radius of convergence greater that the closest packed density p cp . 

We plot the ratios for D > 5 in Figures [21 and In this case the first few virial coefficients are 
positive, and then alternate in sign to the order calculated. We propose the scenario that there is a 
singularity on the positive real axis that dominates the series initially, but at higher order another 
singularity (or singularities) in the complex plane or negative real axis competes with the original 
singularity and hence the new singularity must be at a smaller radius. If the leading singularity 
is on the negative real axis then the ratio plot will smoothly converge to some negative value, 
otherwise the ratios will oscillate. 

For D = 4 in Fig. 01 it seems that despite the absence of negative virial coefficients and the poor 
accuracy of -Bio, an oscillation is developing in the ratio plot in exactly the same way as for D > 5. 
Extrapolation of the series via the methods of Dlog Pades and differential approximants as 
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explained by Guttmann [221, suggests that negative coefficients for D = 4 may occur for k not 
much greater than 12. 

The case D = 3 is plotted in Figures ^ and These ratios do not show the large oscillations 
of D = 4 but close inspection reveals that the slopes are not increasing monotonically as they 
were for D = 2. This may indicate that the plot for D — 3 is displaying very small amplitude 
oscillations, which will eventually result in oscillations in the sign of the coefficients. 




4 Differential Approximants 

We have analyzed the virial coefficients of Table ^ by use of differential approximates using the 
fortran program NEWGRQD given in Guttmann [22]. Our results for the leading singularity on 
the real positive azis in dimensions D — 2,3,4 are tabulated in Table More detailed analysis 
will appear in The notation L,M;N refers to an inhomogeneous first order differential 

approximant, which is the solution of 

zP M {z)f\z) + Q L (z)f(z) = R N (z) (4) 

where the subscript denotes the order of the polynomial, and f(z) is the function that is to be 
approximated. 

One can see from Tables EE1 that there appears to be a singularity on the positive real axis 
close to the space filling density rj = 1 for dimensions D = 2,3,4. The kind of singularity is not 
so clear, for D = 2 there seems to be an algebraic singularity with exponent (j> ~ —1.75, but for 
D = 3, 4 it is not possible to give a good estimate for the exponent. 
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D=s' 
D=6 



0.5 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 

1/k 



Figure 2: Ratio plot for virial coefficients in dimensions D = 5, 6 




0.5 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 

1/k 



Figure 3: Ratio plot for virial coefficients in dimensions D = 7, 8 

5 Conclusion 

In Table ^ above we have reported the first computations of the virial coefficients Bg and Bio for 
hard spheres in dimensions D = 2, ■ • • , 8, and shown that Bg is negative for D > 5. The coefficient 
Bio is negative for D > 4 and at D = 4 the ratios of sucessive coefficients oscillate in such a way 
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3 



2.5 



1.5 





D=4 i — i — i 




/" 







0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 

1/k 

Figure 4: Ratio plot for virial coefficients in dimension D = 4 




0.9 1 1 1 1 1 1 1 1 

0.2 0.18 0.16 0.14 0.12 0.1 

1/k 

Figure 5: Ratio plot for virial coefficients in dimension D = 3 over small domain to show 
non-monotonic behavior of the second derivative. 

as to suggest that negative values may occur for k = 12 — 14. For D = 2, 3 analysis of the first 10 
virial coefficients leads to a radius of convergence greater than close packing. This is in agreement 
with conclusions reached in previous studies based on 8 or fewer virial coefficients. The meaning 
of this is controversial and the present authors argue that the true large k behaviour is not seen 
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Table 3: Singularities from differential approximants on the positive real axis for D = 2,3, 4. 
Blank entries are due to defective approximants. 



Differential 




D = 


z 






9 
O 


LJ — 


A 

~h 


Approximant 








£>2Psing 


<P 


E>2Psing 


<P 


3,3 





1 


987 


-1 


790 


4.068 


-2.818 


7.995 


-3.520 


3,4 





1 


984 


-1 


774 


3.830 


-2.329 


6.843 


-2.478 


4,3 





1 


984 


-1 


774 


3.714 


-2.043 


5.551 


-1.207 


4,4 





1 


987 


-1 


788 


3.732 


-2.090 


7.249 


-2.871 


4,5 





1 


988 


-1 


795 


3.675 


-1.899 


6.985 


-2.583 


5,4 





1 


988 


-1 


795 


3.720 


-2.056 


6.721 


-2.229 


2,2 


1 


1 


946 


-1 


575 


3.659 


-2.014 


6.888 


-2.593 


2,3 


1 

1 


1 


970 


-1 


695 


3.787 


-2.246 


5.412 


-2.922 


3,2 


1 


1 


966 


-1 


677 


4.038 


-2.953 


8.332 


-4.368 


3,3 


1 


2 


021 


-2 


076 


3.811 


-2.298 


7.296 


-2.920 


3,4 


1 


1 


981 


-1 


756 


3.786 


-2.241 


6.764 


-2.389 


4,3 


1 


1 


978 


-1 


740 


3.708 


-2.024 


4.663 


-3.708 


2,1 


2 


1 


945 


-1 


572 


3.676 


-2.054 






2,2 


2 


1 


967 


-1 


682 


3.641 


-1.987 


7.715 


-3.531 


2,3 


2 


2 


008 


-1 


900 


3.799 


-2.269 


7.277 


-2.916 


3,2 


2 










3.874 


-2.468 


7.250 


-2.899 


3,3 


2 


1 


971 


-1 


679 


3.773 


-2.210 


7.308 


-2.920 


2,1 


3 


1 


959 


-1 


628 


3.599 


-1.889 


6.802 


-2.680 


2,2 


3 


1 


984 


-1 


784 


3.747 


-1.956 


7.508 


-3.404 


2,3 


3 


1 


982 


-1 


778 


3.777 


-2.110 


7.003 


-2.539 


3,2 


3 


1 


981 


-1 


770 


3.779 


-2.034 







Table 4: Approximate position of singularities with exponents 



D 


B2Psing 






1 


1.00 


1.00 


-1.00 


2 


1.98 


0.99 


-1.75 


3 


3.75 


0.94 


-2.1 


4 


7.00 


0.88 


-3 



in the first 10 coefficients. More complete analysis and discussion of this is given in |19j . 
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